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Large  crystalline  molecular  shells  buckle  spontaneously  into  icosahedra  while 
multicomponent  shells  buckle  into  various  polyhedra.  Continuum  elastic  theory 
explains  the  buckling  of  closed  shells  with  one  elastic  component  into  icosahe¬ 
dra.  A  generalized  elastic  model,  on  the  other  hand,  describes  the  spontaneous 
buckling  of  inhomogeneous  shells  into  regular  and  irregular  polyhedra.  By  co¬ 
assembling  water-insoluble  anionic  (-1)  amphiphiles  with  cationic  (3+)  amphiphiles, 
we  realized  ionic  vesicles.  Results  revealed  that  surface  crystalline  domains  and 
the  unusual  shell  shapes  observed  arise  from  the  competition  of  ionic  correla¬ 
tions  with  charge-regulation.  We  explain  here  the  mechanism  by  which  these  ionic 
membranes  generate  a  mechanically  heterogeneous  vesicle.  ©  2016  Author(s).  All 
article  content,  except  where  otherwise  noted,  is  licensed  under  a  Creative 
Commons  Attribution  (CC  BY)  license  (http://creativecommons.Org/licenses/by/4.0/). 
[http://dx.doi.org/10. 1063/1.4953570] 


Convex  polyhedral  shells  have  been  observed  in  various  microscopic  systems  including  viral 
capsids,1  fullerenes,2  protein-based  bacterial  organelles,  '  and  co-assembled  ionic  amphiphiles.4  The 
most  commonly  observed  polyhedral  symmetry  in  self-organized  homogeneous  structures  is  the 
icosahedron,  which  has  the  highest  possible  symmetry.  We  explain  here  a  way  to  design  various 
polyhedral  shapes  with  reduced  symmetry  using  elasticity  theory5  and  full  atom  and  coarse-grained 
computer  simulations  with  input  from  experimental  studies.4 

Elasticity  theory6  explains  that  the  lowest-energy  configuration  of  a  buckled  crystalline  shell 
is  an  icosahedron.7  The  theory  describes  how  any  spherical  crystalline  shell  made  of  a  homoge¬ 
nous,  isotropic,  elastic  material  buckles  into  an  icosahedron  when  the  Foppl-von  Karman  number 
y  =  YR2/ k,  where  Y  is  the  Young’s  modulus,  k  is  the  bending  rigidity,  and  R  is  the  linear  size  of  the 
shell,  is  sufficiently  high.  Since  closed  crystalline  membranes  have  at  least  12  five-fold  disclinations 
that  repeal  each  other,  these  defects  are  positioned  at  the  vertices  of  an  inscribed  icosahedron  in 
the  sphere.8  When  the  Foppl-von  Karman  number  exceeds  a  critical  value,  these  defects  buckle 
simultaneously,  and  the  spherical  shell  acquires  an  icosahedral  shape. 

Recently,  we  found  that  elastic  membranes  with  elastic  heterogeneities,  which  are  achieved 
via  increasing  the  number  of  the  shell  components,  buckle  into  regular  and  irregular  polyhedra, 
as  shown  in  Figure  l.5  We  analyzed  the  simplest  case  of  a  two-dimensional  membrane  composed 
of  two  components  A  and  B  that  are  individually  isotropic  and  characterized  by  their  Young’s 
modulus,  Ya ,  Yb,  and  bending  rigidity  ka,  kb;  the  Poisson  ratios  were  set  equal  to  vA  =  vB  ~  0.33. 
We  found  the  lowest  energy-state  configuration  by  discretizing  a  model  for  elastic  membranes.9  The 
Hamiltonian  is  given  by  the  sum  of  two  contributions,  the  bending  energy  Eb,  and  the  elastic  Es 
energy  given.  The  bending  energy  is 
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FIG.  1.  Snapshots  of  configurations  at  different  fraction  f.  The  A-component  (in  blue)  has  YA  =  5.77  and  k'A  =  0.06  while 
the  B-component  (in  red)  has  Yg  =  5.77  and  Kg  =  28.9.  The  radius  of  the  shell  is  about  R  ~  1 1.3  in  units  of  the  average  bond 
length.  The  fraction  f  varies  (from  left  to  right,  top  to  bottom)  as  f  =  0%,  4%,  10%,  20%,  30%,  40%,  50%,  60%,  65%,  70%, 
75%,  80%,  85%,  90%,  95%,  95.5%,  96%,  96.5%,  97%,  97.5%,  97.7%,  98%,  98.5%,  99%,  100%.  Adapted  with  permission 
from  Vemizzi  eta/.,  Proc.  Natl.  Acad.  Sci.  U.  S.  A.  108(11),  4292  (2011).  Copyright  2011  National  Academy  of  Sciences, 
USA. 


Eb  =  J  dS  [2 K(r)(H  -  H0 f  +  kq{t )Gj  , 

where  H  is  the  mean  curvature,  and  G  is  the  Gaussian  curvature,  and  the  bending  (k)  and  Gauss  (kg) 
rigidities  are  a  function  of  the  position  in  the  shell  surface  r(si,S2).  The  elastic  energy  is 

Es=l-JdS  [,l(Tr(w))2  +  2/iTr  (u2)]  , 

where  X,  /;  are  the  position  dependent  Lame  coefficients  (for  each  component,  Xl-/ji  =  Y|3/8, 
where  Y  is  the  Young  modulus  of  component  I  =  A,  B).  The  lowest  energy  configuration  was 
obtained  by  Monte  Carlo  simulated  annealing  simulations  for  a  fixed  relative  fraction  /  of  the  two 
components. 

Possible  buckled  structures  can  be  obtained  considering  components  A ,  B  with  intermediate 
compositions  /  of  the  B  component  such  that  the  stable  shape  in  the  limits  of  homogeneous  compo¬ 
sition  is  icosahedral  (/  =  0)  or  spherical  (/  =  1),  respectively.  A  sufficient  condition  for  that  to 
happen  is  to  choose  Kb  >  Ka  such  that  Yb  <  y*  <  Y a>  where  Ya  =  YaR2/ka  and  ys  =  YbR2/kb-  We 
considered  the  B-component  more  rigid  than  the  A-component  (kb  >  Ka). 

At  intermediate  values  of  the  composition  /,  the  component  with  low-bending  rigidity  tends  to 
form  lines  that  merge  into  vertices.  The  pattern  of  lines  aids  the  faceting  of  the  shell  into  a  polyhedral 
structure  often  with  three-fold  vertices.  Such  vertices  are  joined  together  by  sharp  edges  occupied 
by  component  A  that  surrounds  polygons.  Even  though  the  components  are  chemically  compatible. 
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they  segregate  due  to  their  different  elastic  properties.  This  can  be  rationalized  using  geometric  argu¬ 
ments.  1,1  The  Gauss-Bonnet  theorem  states  that  the  integral  of  the  Gaussian  curvature  is 


where  D  is  the  integral  over  the  domain,  <9D  is  the  domain  boundary,  kg  is  the  geodesic  curva¬ 
ture,  and  Xd  is  the  Euler  characteristic  of  the  surface.  Therefore,  applying  this  theorem  to  the 
non-homogenous  closed  shell  with  soft  and  hard  domains,  one  obtains  that  the  Gaussian  contribu¬ 
tion  to  the  energy  is,10 


where  we  used  the  identity  that  equates  the  sum  of  the  Euler  characteristics  of  the  soft  and  hard 
domains  to  the  Euler  characteristic  of  a  closed  surface,  which  is  two  (/soft  +  Xhard  =  Xd  =  2)  and 
also  that  <9D|lal-d  =  -<9Dsoft,  since  the  boundary  of  the  soft  component  is  the  boundary  of  the  hard 
component  but  in  the  opposite  direction.  Given  that  in  the  discretized  model9  used  to  generate 
Figure  1,  the  Gaussian  rigidity  and  the  bending  rigidity  are  related  by  Kg1  =  -2k1  for  both  hard  and 
soft  components  (I  =  soft,  hard)  we  obtain, 


Since  Khard  >  Ksof’,  both  terms  favor  a  large  negative  <9Dsoft,  and  a  large  positive  integral  geodesic 
curvature  of  the  boundary.  This  means  that  a  soft-component  domain  tends  to  have  several  bound¬ 
aries,  Xsoft  =  2  —  h,  where  h  is  the  number  of  boundaries  for  a  connected  domain,  and  to  be  negative, 
h  should  be  at  least  three.  We  note  that  the  Euler  characteristic  of  a  disk  is  one.  Consequently,  the 
hard-component  domain  tends  to  fragment  in  several  disconnected  domains,  roughly  h  in  number, 
each  one  with  the  topology  of  a  disk. 

The  Gauss-Bonnet  theorem  stresses  the  numerical  results  that  show  the  buckling  in  the  case 
of  mechanically  heterogeneous  shells  is  not  mediated  by  defects,  as  in  the  case  of  homogenous 
crystalline  shells.  That  is,  even  if  the  Hamiltonian  for  the  heterogeneous  shells  is  discretized  using 
a  triangular  mesh  with  fixed  connectivity  and  with  the  required  12  five-fold  disclinations,  those 
defects  are  located  in  the  flat  regions  of  the  shell  (and  are  slightly  sunken  in).  The  buckling  into 
polyhedra  occurs  by  segregation  of  the  component  with  lowest  bending  rigidity  into  narrow  bands 
that  fold;  that  is,  by  concentrating  the  curvature,  which  is  a  mechanism  that  also  describes  the 
crumpling  of  membranes.11 

We  note  that  in  the  numerical  model,  polyhedral  structures  are  expected  only  in  specific  regions 
of  the  parameter  space  of  the  two-component  system  (R,  Y, t,  k aYb,  k b)-  Shape  competition  between 
the  elastic  components  is  required  for  buckling  into  polyhedral  shells.  Faceting  in  a  shell  with  a 
given  radius  occurs  when  at  least  one  of  the  components  has  a  large  bending  rigidity  and  a  low 
Young’s  modulus.  For  very  small  shells,  however,  an  exceedingly  high  Young’s  modulus  and  a  very 
low  bending  rigidity  are  required.  Therefore,  faceting  effects  are  expected  to  occur  for  a  narrow 
range  of  shell  sizes  where  the  shape-competition  is  consistent  with  the  elastic  parameters  of  soft 
materials  (at  least  one  component  would  have  to  have  either  extremely  large  or  extremely  low 
Y-to-K  ratio). 

Regardless  of  the  apparently  highly  narrow  range  of  parameter  space  required  for  the  reali¬ 
zation  of  the  buckling  of  mechanically  heterogeneous  elastic  shells  into  polyhedra,  the  model  is 
very  robust.  As  discussed  above,  the  Gauss-Bonnet  theorem  supports  fragmentation  of  the  shells 
driven  by  the  formation  of  ridges  by  the  component  of  lowest  bending  rigidity.  The  formation  of 
ridges  in  the  buckling  of  membranes  has  long  been  established. 12  Therefore,  it  is  not  surprising  that 
such  effect  escapes  closed  membranes’  buckling  when  the  folding  is  facilitated  (i.e.,  by  adding  a 
second  component  that  can  segregate  into  ridges).  Moreover,  it  has  been  shown  that  faceted  bilayer 
vesicles  that  allow  the  segregation  of  excess  amphiphiles  along  the  ridges  of  polyhedra  do  have 
lower  bending  energies  than  spherical  bilayer  vesicles.13  Therefore,  the  required  ingredients  for 
buckling  into  polyhedral  shells  are  realizable.  Indeed,  experiments  on  bilayers  of  co-assembled 
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cationic  and  anionic  amphiphiles  revealed  giant  icosahedral  vesicles.14  These  giant  icosahedra 
have  pores  at  the  12  vertices  of  the  faceted  vesicles  and  are  only  formed  at  zero  salt  conditions 
when  there  is  an  excess  of  an  amphiphile  component,  which  segregates  at  the  opening  of  the 
pores. 

We  recently  demonstrated  the  buckling  into  various  polyhedra  without  pores  by  co-assembling 
amphiphiles  with  various  degrees  of  ionization  into  ionic  vesicles.4  We  generated  crystalline  mem¬ 
branes  in  a  large  variety  of  geometries,  including  faceted  polyhedral  shells  by  co-assembling 
cationic  (+3)  and  anionic  (-1)  amphiphiles  at  different  pH  values.  Although  single-tailed  amphi¬ 
philes  typically  form  micelles  in  aqueous  solutions,  co-assembled  cationic-anionic  mixtures  can 
form  bilayers.15  In  fact,  work  with  polymerizable  amphiphiles16  showed  that  a  large  charge  imbal¬ 
ance  between  the  cationic  and  anionic  head  groups  enables  their  co-assembly  into  nanoscale 
buckled  vesicles.  By  co-assembling  oppositely  charged  surfactants  without  a  polymerizable  group, 
we  demonstrated  that  electrostatics  drives  the  crystallization  of  the  tails  and  determines  the  crys¬ 
talline  structure  of  the  domains  resulting  into  buckled  nanoscale  vesicles  that  are  stable  at  high 
salt  concentrations.4  The  crystalline  domains,  controlled  by  pH,  lead  to  the  formation  of  shapes 
other  than  icosahedra,  very  similar  to  the  shapes  obtained  in  the  polymerized  co-assembled  vesi¬ 
cles.16  These  shapes  resemble  organelle  microcompartment  geometries3  and  the  faceted  shape  of 
halophilic  organisms.17 

The  selected  3+  cation  C16-K3  dissolves  in  water  and  forms  micelles  of  ~10  nm  diameter. 
However,  it  can  be  co-assembled  with  water-insoluble  palmitic  acid  C16-COOH  at  different  pH 
values.  The  TEM  images  of  the  cation/anion  mixture  of  amphiphiles  show  a  pH  dependent  diversity 
of  shapes  including  faceted  vesicles  and  ribbons.4  Collected  in  situ  small  angle  x-ray  (SAXS)  and 
wide  angle  x-ray  (WAXS)  data  demonstrated  that  the  structures  were  consistent  with  a  bilayer 
crystalline  membrane  morphology.  The  WAXS  data  showed  for  pH  <  7  and  pH  >10,  the  appear¬ 
ance  of  a  single  diffraction  peak,  which  indicates  that  the  molecules  are  packed  in  a  2D  hexagonal 
lattice;  whereas  the  presence  of  two  diffraction  peaks  for  pH  8-9  indicates  a  rectangular-C  structure. 
Vesicle  shapes  were  observed  by  TEM  at  pH  5  and  pH  10  while  at  pH  8,  nanoribbons  and  sheets 
were  observed.  These  experimental  results  demonstrated  that  changes  on  the  mesoscale  structure 
are  due  to  the  nanoscale  packing  of  the  molecules. 

Monte  Carlo  simulations  were  used  to  determine  how  pH  relates  to  the  degree  of  ionization 
of  the  molecules.18  Here  we  summarize  the  main  results  of  the  Monte  Carlo  simulations  of  a 
charge-fluctuating  model  that  treats  charge  correlations  in  detail.  The  degree  of  dissociation  of  a 
lipid  bilayer  depends  on  both  the  salt  concentration  and  the  pH  value.  If  the  surface  is  acidic,  that 
is  a  dissociated  head  group  has  a  negative  charge,  at  low  pH,  the  excess  H+  ions  in  the  solution 
will  inhibit  the  dissociation  and  the  bilayer  will  be  essentially  neutral.19  As  the  pH  increases,  the 
dissociation  becomes  more  favorable  and  the  surface  gradually  becomes  charged.  When  pH  is 
equal  to  the  pKa  of  a  given  molecule,  on  average  50%  of  the  head  groups  are  dissociated.  When 
interactions  between  charges  on  the  bilayer  surface  are  included,  the  electrostatic  repulsion  between 
two  negative  charges  will  inhibit  the  dissociation.  Meanwhile,  when  a  positive  charge  is  next  to  the 
negative  charge,  the  average  dissociation  rate  will  be  enhanced. 

We  modeled  the  effects  of  electrostatic  interactions  via  a  discrete,  lattice  model  introduced  by 
Netz.18  We  constructed  an  appropriate  lattice  with  lattice  constant  a  =  1,  that  is,  use  the  lattice  spac¬ 
ing  as  the  unit  of  length.  In  our  particular  experimental  system  of  interest,  a  roughly  corresponds  to 
0.5  nm.  The  Hamiltonian  is  given  by 


H 

knf 


A 'A 

kBT 


I<+02».s+52>1>",i 


iK 


where  77,7  =  0  (1)  if  sites  i  and  j  are  of  the  same  (different)  type  and  nf(B)  -  0  or  I  is  the  dissocia¬ 
tion  state  of  the  acidic  (basic)  site  i.  The  Debye  screening  length  is  given  by  £  =  1/(4  nlB  2z?C;)^2, 
where  n  and  c,-  and  valence  and  the  concentration  of  salt  ions,  and  the  Bjerrum  length,  /  B  = 
ell‘Xns{)ErkHT ,  is  the  length  at  which  electrostatic  interaction  is  comparable  with  thermal  energy 
k bT,  where  e  is  the  unit  charge,  sr  is  the  relative  permittivity  of  the  medium,  and  eq  is  the  dielectric 
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permittivity  of  vacuum.  The  chemical  potentials  are  given  by 

/-rp  =  -  In  10  (pH  -  pKci(A)) 

~  In  10  (pH  -  pKa(B)) 

Monte  Carlo  simulations  were  used  to  determine  the  average  dissociation  fraction  of  the  co¬ 
assembled  system  as  a  function  of  pH  and  Debye  screening  length  £.  While  the  charge  dissociation 
for  single-component  systems  is  heavily  suppressed  by  the  electrostatic  repulsion  between  neigh¬ 
boring  sites  (ideal  case),  in  the  case  of  mixtures,  the  dissociation  is  favorable  and  the  system  is 
completely  charged  over  a  range  of  pH  values. 

Molecular  dynamics  (MD)  simulations  were  performed  at  a  low  degree  of  ionization  of  the 
palmitic  acid  groups  to  mimic  the  low  pH  data  and  a  high  degree  of  ionization  to  explain  the 
high  pH  data.  The  atomistic  simulations  were  carried  out  for  electroneutral  systems  with  30% 
(Figure  2(a))  and  95%  (Figure  2(b))  average  ionization  of  the  palmitic  acid  molecules,  which  corre¬ 
spond  roughly  to  pH  4  and  8,  respectively.  Interdigitation  of  the  upper  and  the  lower  layer  in  the 
bilayer  was  observed  only  in  the  system  with  95%  ionization  of  palmitic  acid.  The  rectangular-C 
packing  of  the  tails  observed  experimentally  was  not  reproduced  by  the  MD  simulations.  This  is 
possibly  due  to  the  force  fields,  which  were  not  parameterized  for  crystalline  phases. 

A  simple  estimate  of  the  electrostatic  correlations  can  explain  the  formation  of  the  crystal¬ 
line  phase  in  the  vesicles.  The  head  groups  located  between  polar  water  (er  ~  80)  and  nonpolar 
tails  (sr  as  1)  experience  a  dielectric  response  of  the  order  of  the  mean  between  the  two  me¬ 
dia,  er  ~  40,  resulting  in  1r  e2/4nsosrkBT  ~  1 .4  nm.  The  electrostatic  energy  to  hold  together 


FIG.  2.  Cationic  and  anionic  amphiphiles  design  to  co-assemble  and  organize  into  polyhedral  shell  shapes  determined  by 
a  combination  of  WAXS-SAXS,  TEM  (not  shown),  full  atom  and  coarse  grained  molecular  dynamic  (MD)  simulations, 
and  continuum  elasticity.  Bilayers  at  pH  4  (a)  and  at  pH  8  (b)  from  full  atom  MD  using  the  GROMACS  4.5.5  package, 
(c)  Vesicles  with  shapes  expected  for  heterogeneous  elastic  closed  membranes  at  low  pH  values  from  coarse  grained  MD 
simulations  using  the  MARTINI  version  2.1  force  field,  (d)  The  continuum  elastic  model  of  mechanically  heterogeneous 
shells  (hard  and  soft  components  are  shown  in  yellow  and  blue,  respectively)  shows  that  for  fixed  crystalline  domain  size, 
small  vesicles  (left:  Khard/Ksoft  =  50  for  95%  of  hard  and  5%  of  components)  have  low  symmetry  shapes  while  large  vesicles 
(right:  Khard/KSoft  =  20  for  85%  of  hard  and  15%  of  soft  components)  have  more  spherical  shapes  in  agreement  with  TEM 
images.4  Insets:  vesicle  height  map  with  regions  shown  in  red  (blue)  is  closest  (furthest)  to  (from)  the  center  of  the  vesicle; 
the  ratio  of  the  largest  (Rmax)  to  the  lowest  (Rmin)  radius  is  Rmax/Rmin~  1.1.  WAXS  shows  that  the  alkyl  tails  packing 
changes  with  pH  from  hexagonal  (pH  4-7)  to  rectangular-C  (pH  8-9.5)  to  hexagonal  lattices  (pH  10  and  above).  Adapted 
with  permission  from  Leung  et  al.,  ACS  Nano  6(12),  10901  (2012).  Copyright  2012  American  Chemical  Society. 
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one  +3  and  three  -1  chains  into  a  four-tail  bundle  is  then  Ee  ~  3(-l)lBkBT/d  ~  -9kBT,  where 
we  have  used  the  tail  inter-chain  distance  observed  by  WAXS  ( d  ~  0.5  nm).  This  estimate  is 
in  agreement  with  the  estimates  from  atomistic  MD  simulations  of  about  -2.8  kBT  (30%)  per 
NH3+  -COO-  ion  pair  for  their  cohesive  energy.4  Adding  the  hydrophobic  attraction  from  atom¬ 
istic  MD  simulations, -20A:Br/chain,  gives  A’/,  ~  -90 kBT.  In  contrast  to  ionic  driven  assembly  of 
incompatible  molecules20  and  of  molecules  without  reinforced  attractions,21  the  ionic  driven  assem¬ 
bly  discussed  here  induced  short-range  attractions  among  dissimilar  molecules  of  opposite  charge, 
explaining  the  observed  stability  of  the  ionic  crystalline  bilayer  with  the  addition  of  salt  up  to  salt 
concentration  n  =  500  mM  NaCl. 

In  order  to  reveal  the  structure  of  the  faceted  vesicles  at  low  pH,  simulations  with  a  MARTINI 
coarse-grained  MD  model22  were  also  performed.  In  Figure  2(c),  a  MARTINI  MD  simulation  snap¬ 
shot  shows  that  the  internal  structure  of  the  bilayer  in  the  vesicle  is  not  homogenous.  The  vesicle 
can  be  thought  of  as  a  two-component  shell,  where  one  component  represents  the  crystalline  do¬ 
mains  at  the  vesicle  faces  and  the  second  component  represents  the  liquid-like  boundaries  between 
the  hard  facets  where  the  curvature  is  large.4  The  bilayer  packing  and  thickness  in  the  flat  areas  are 
in  agreement  with  the  full  atom  MD  simulations  (see  Figure  2(a)). 

Given  that  the  size  of  the  crystalline  domains  obtained  by  SAXS  is  about  30  x  30  nm,2  the 
small  vesicles  of  100  nm  to  200  nm  in  diameter  have  only  few  crystalline  domains  covering 
their  surfaces.  Therefore,  these  are  more  likely  to  be  faceted  into  polyhedral  structures  with  low 
symmetry;  meanwhile  large  vesicles  have  nearly  spherical  faceted  shapes.  Both  the  small  highly 
asymmetric  and  the  large  nearly  spherical  vesicles  that  were  observed  by  TEM  were  recovered  by 
the  two-component  elastic  model  (see  Figure  2(d)).  The  crystalline  domains  of  the  vesicle  faces 
represent  the  hard  component  with  a  bending  modulus  Khard,  and  the  second  (soft)  component 
models  boundaries  between  the  hard  facets,  which  are  thinner  and  softer  with  a  bending  modulus 
Ksoft.  The  shapes  in  Figure  2(d)  show  vesicle  shapes  obtained  by  simulated  annealing  Monte  Carlo 
optimization  of  a  discrete  triangulation4  for  large  and  small  crystalline  domains  sizes,  respectively 
(Yimrd/Ysoft  =  1  since  the  Young  modulus  is  less  sensitive  to  thickness  variations  than  the  bending 
rigidity2  ’  and  a  related  model  of  mechanically  heterogeneous  elastic  shells24  with  Yharci/Ysof,  >  1 
gives  similar  results).  We  note  that  if  the  edges  of  the  polyhedra  are  liquid-like,24  the  shapes  are 
not  likely  to  change  if  the  fraction  of  fluid  to  crystalline  phases  is  low.  Interestingly,  the  phase 
diagram  for  tail  packing,  as  observed  in  Langmuir  films, 24-2h  allows  fluid-hexagonal  coexistence. 
However,  there  does  not  exist  an  obvious  coexistence  region  between  rectangular-C  and  fluid  phases 
in  monolayers,  which  may  explain  why  the  vesicles  are  only  observed  at  the  pH  values  when  the 
tails  are  crystallized  in  a  hexagonal  lattice  and  are  not  stable  at  pH  is  8  when  the  tails  crystalline 
structure  is  rectangular-C. 

Interestingly,  the  MARTINI  MD  simulations  show  disordered  regions  with  high  mean  curva¬ 
tures  as  predicted  by  the  continuum  model  of  heterogeneous  elastic  shells,  even  though  these 
liquid-like  regions  have  very  low  or  null  young  moduli  (Fig.  3).  However,  since  the  Gauss-Bonnet 
analysis  only  requires  large  bending  rigidity  differences  (not  Young  modulus  differences),  the  frag¬ 
mentation  of  shells  should  be  observed  in  shells  with  liquid  and  solid  domains,  provided  the  fraction 
of  surface  area  covered  by  the  liquid-like  regions  is  very  small  compared  to  the  fraction  of  area 
covered  by  the  crystalline  regions  such  that  these  easy  to  bend  regions  can  form  ridges. 

In  conclusion,  the  low  symmetry  faceted  shapes  in  co-assembled  vesicles  of  oppositely  charged 
single-tail  amphiphiles  arise  due  to  crystalline  domains  induced  by  ionic  correlations  among  the 
headgroups  and  van  der  Waals  interactions  between  the  tails.  These  ionic  correlations,  as  in  ionic 
bonding,  are  not  screened  even  at  high  salt  concentrations  generating  membranes  stable  up  to 
500  mM  of  NaCl.  Their  crystalline  lattice  and  morphology,  however,  are  determined  by  the  solution 
pH.  The  observed  irregular  faceted  shapes  arise  only  in  hexagonal  crystalline  lattices,  where  crys¬ 
talline  domains  are  separated  by  soft  interfaces  that  bend  to  release  stress.  Since  ionic  crystalline 
domains  also  arise  upon  the  adsorption  of  multivalent  ions  on  lipids,27  our  results  suggest  that 
spherical-to-faceted  reversible  transitions  can  be  induced  in  cationic  lipid  vesicles  by  adding  (or 
chelating)  metallic  multivalent  ions. 

In  the  present  study,  we  found  that  for  crystalline  domains  of  the  same  sizes,  small  vesicles  are 
more  likely  to  be  faceted  into  polyhedral  structures  with  low  symmetry  (lower  than  the  icosahedra) 
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The  mean 
permittivity 

4=14  A 


For  3:1  with  d~  5A 

ec(d)/KBT~-7AI/d~9 


Hydration  of  4  tails  ~  80  KBT 
Cohesive  energy  of  a  neutral  3:1  unit  ~90  KBT 


FIG.  3.  Schematic  of  the  cross-sectional  area  and  top  view  of  the  co-assembled  amphiphiles  configuration  that  leads  to  the 
tail  crystallization  using  the  parameters  from  the  x-ray  experiments  and  estimates  from  the  full  atom  simulations;  here,  ec  is 
the  electrostatic  energy  of  a  neutral  3:1  unit. 


and  large  vesicles  in  which  defects  spread  over  many  boundaries  have  nearly  spherical  faceted 
shapes.  The  kinetic  path  for  the  assembly  of  these  faceted  polyhedral  vesicles  may  resemble  the 
initial  kinetic  path  observed  in  the  formation  of  liquid  vesicle  where  disk-like  intermediate  micelles 
are  observed. 2S  In  the  electrostatic  driven  co-assembly  case,  however,  crystalline  bilayers  of  various 
shapes  may  develop  initially  that  subsequently  merge  into  closed  structures. 

This  combined  theoretical,  numerical,  and  experimental  work  explains  the  relation  of  the  crys¬ 
talline  structure  of  the  membrane  to  their  geometries  and  elucidates  the  role  of  ionic  interactions  in 
the  assemblies.  The  closed  shells  are  shown  to  be  stable  only  when  the  membrane  crystal  symmetry 
is  hexagonal.  This  discovery  will  help  researches  to  design  shells  in  various  geometries  and  eluci¬ 
dates  the  mechanisms  used  in  nature  to  self-assemble  cellular  microcompartments  and  the  shapes 
of  halophilic  organisms.  For  example,  the  robust  envelope  of  Walby’s  rectangular  bacterium  takes 
a  square  geometry  in  natural  brines,  even  though  the  cell  wall  is  actually  composed  of  a  hexagonal 
lattice29  of  various  lipids. 

In  summary,  we  have  described  the  mechanism  that  connects  the  mesoscale  structure  of  a  com¬ 
plex  system  to  the  molecular  scale  conformation  of  the  components.  We  exploited  such  connections 
to  design  a  family  of  polyhedral  shells.  The  elements  required  in  the  design  involved  geometric 
principles,  elasticity  theory,  computer  simulations,  and  experimental  results.  This  example  demon¬ 
strates  that  robust  multicomponent  structures  can  be  assembled  when  ionic  correlations  reinforce 
molecular  interactions  at  the  nanoscale  such  as  the  packing  of  hydrophobic  tails  into  crystalline 
lattice.  In  this  case,  the  ionic  assemblies  are  stable  at  high  salt  concentrations  since  the  assemblies 
are  self-screened.  The  mesoscale  structure  is  a  direct  consequence  of  the  constraints  imposed  by  the 
geometrical  packing. 
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